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RECONSTRUCTION OF ENERGY-DEPENDENT 
STURM-LIOUVILLE EQUATIONS FROM TWO 

SPECTRA 

NATALIYA PRONSKA 



>>'. 

d , Abstract. In this paper we study the inverse spectral problem of 

reconstructing energy-dependent Sturm-Liouville equations from 
two spectra. We give a reconstruction algorithm and establish ex- 
istence and uniqueness of reconstruction. Our approach essentially 
exploits the connection between the spectral problems under study 
and those for Dirac operators of a special form. 



Oh 

1. Introduction 

In this paper, we consider energy-dependent Sturm-Liouville differ- 
ential equations 

(1.1) - y" + qy + 2Xpy = \ 2 y 

G\ ; on (0, 1); here A G C is the spectral parameter, p is a real- valued func- 

tion from £2(0, 1) and q is a real- valued distribution from the Sobolev 
space W / 2~ 1 (0, 1), i.e. q = r' with a real- valued r E L 2 (0, 1) (see a de- 
tailed definition in the next section). Our aim is to study the inverse 
CN ' problem of reconstructing the potential p and a primitive r of q from the 

spectra of the problem (II -ip under two types of boundary conditions: 
the Dirichlet ones 

(1.2) v(0)=y(l) = 

en ' 

and the so-called mixed conditions 

(1.3) y (0) = (y'-ry)(l) = 0. 

The spectral problems under consideration often arise in classical 
and quantum mechanics. The most common example is modelling 
of the motion of massless particles such as photons by means of the 
Klein-Gordon equations, which can be transformed to Sturm-Liouville 
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equations with potentials depending on the spectral parameter. The 
corresponding evolution equations are used to describe interactions of 
colliding spinless particles. The equations of interest also arise in mod- 
elling of mechanical systems vibrations in viscous media, see [16]. As 
the spectral parameter is related to the energy of the system, this ex- 
plains the terminology "energy-dependent potentials" widely accepted 
in the physical and mathematical literature. 

Since equation (II. ip depends nonlinearly on A, we should regard the 
spectral problems (II. ip . (II. 2p and (II. ip . (II. 3p as those for quadratic 
operator pencils. Some spectral properties of such problems can be 
derived from the general spectral theory of polynomial operator pen- 
cils [31]. 

An interesting approach to the spectral analysis of Klein-Gordon 
equations, using the theory of Krein spaces (i.e. spaces with indefinite 
scalar products) was suggested by P. Jonas [22] and H. Langer, B. Na- 
jman, and C. Tretter [26l[27J[36] . This method was also applied to the 
spectral analysis of problems of interest, cf. [58] . 

Equation (ll.ip has appeared earlier in the context of inverse scat- 
tering. For example, M. Jaulent and C. Jean studied inverse scat- 
tering problems for energy-dependent Schrodinger operators on the 
line and half-line in [T8Tl2~T] . Such problems were also considered 
in [Il[23l[29l[30l[32l[33l[39l[lll[l5] . 

However, the inverse problems for energy- dependent Sturm-Liouville 
equations on a finite interval have not been studied sufficiently well. 
M. Gasymov and G. Guseinov discussed such problems with p e 
W^iO, 1) and q G L 2 (0, 1) and with Robin boundary conditions in 
their short paper [SJ of 1981 containing no proofs. In the papers 
[SICEnilMlESJIlI] the spectral problems of interest with (quasi)-periodic 
boundary conditions were considered, but only the Borg-type unique- 
ness results were obtained therein. 

In the present paper, we consider the inverse spectral problem 
for ( II .ip under minimal smoothness assumptions on the real-valued 
potentials p and q, including, e.g. the case when q contains Dirac delta- 
functions and/or Coulumb-like singularities. Such functions often arise 
in modelling of interactions in molecules and atoms in quantum me- 
chanics, see the monographs by S. Albeverio et al. [2J and by S. Albev- 
erio and P. Kurasov [I] and the extensive reference lists therein. The 
inverse spectral problem for equation (II. ip under such minimal smooth- 
ness assumptions on the potentials p and q was considered in [H] , where 
the reconstruction from the Dirichlet spectrum and the set of suitably 
defined norming constants was investigated. 
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Here, we study the problem of reconstruction of the potentials p and 
q of (II. ip from the two spectra determined respectively by the bound- 
ary conditions (11. 2p and (11.31) . Namely, we give a complete description 
of these spectra; establish constructively existence of the pencil with 
given spectra and thus derive a reconstructing algorithm; finally, we 
prove uniqueness of the reconstruction. Note that (I1.3P contains an un- 
known primitive r of q; therefore it is in fact this primitive r (called the 
regularized potential, or simply the potential) rather than just q that 
should be found. Our approach uses the technique analogous to that 
of P3|. Namely, we shall associate the spectral problems (II. ip . (II. 2p 
and ( 11. 1 p . ( 11. 3p with those for Dirac operators of a special form and 
shall essentially use the well-developed inverse spectral theory for Dirac 
operators in our analysis. 

The paper is organized as follows. In the next section we give the ex- 
plicit formulation of the inverse problem of interest and state the main 
results. In Section [31 we show that the spectral problems under con- 
sideration can be reduced to those for Dirac operators with potentials 
of a special form. In Section HJ we introduce the so-called transforma- 
tion operators relating such Dirac operators to the Dirac operators in 
a canonical form. Existence and uniqueness of Dirac operators with 
potential of the special form having given sets as their spectra is estab- 
lished in Section [5j in this section we also prove the main results. In 
the last section we formulate the reconstruction algorithm and discuss 
some possible extensions. 

Notations. Throughout the paper, we denote by L 2 ,ir(0, 1) and 
^"^(0, 1) the spaces of real- valued functions in L 2 (0, 1) and distribu- 
tions in W 2 ~ 1 (0, 1), respectively, and by Ai 2 = A^C) the linear space 
of 2 x 2 matrices with complex entries endowed with the Euclidean 
operator norm. Next, p(T) and o~(T) will stand for the resolvent set 
and the spectrum of a linear operator or a quadratic operator pencil T. 
The superscript t will signify the transposition of vectors and matrices, 
e.g. (ci,c 2 y is the column vector (^). 



2. Preliminaries and main results 

As was mentioned before, the spectral problems (II. ip . (II. 2p and (II. ip . 
( II. 3p are the spectral problems for quadratic operator pencils. To give 
a precise definition, we consider first the differential expression 

t(y) ■= -y" + qy- 
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Recall that q = r' is a real- valued distribution in the Sobolev space W 2 ~ 1 (0, 1), 
and to define £ rigorously one can use, e.g. the regularization by quasi- 
derivative method due to Savchuk and Shkalikov [401I4T] . Following 
that method, one introduces the quasi- derivative y^ := y' — ry of an 
absolutely continuous y and then define t as 

t{y) = -(y m )'-ryW-r*y 

on the set 

domi = {ye AC(0, 1) | E AC[0, 1), i{y) E L 2 (0, 1)}. 

It is easy to verify that £(y) = —y" + qy in the sense of distributions 
and that for regular potentials q the above definition coincides with the 
standard one. 

Denote by A\ and A 2 Sturm-Liouville operators defined via 

A 3 -y = £(y), j = l,2, 

on the domains 

dom^! := {y E domi | y(0) = y [1] (l) = 0} 

and 

domv4 2 := {y E domi \ y(0) = y(l) = 0}. 
For q E W^r(0, 1) the operators A% and A 2 are known |101I1T] to be 
self-adjoint, bounded below, and to have simple discrete spectra. 

Further, we denote by B the operator of multiplication by the po- 
tential p E L 2 (0, 1). The operator B is in general unbounded. How- 
ever, since the Green functions of the operators Aj are continuous on 
the square [0,1] x [0,1], the resolvents of Aj are of Hilbert-Schmidt 
class. It follows that for every non-real A the operators B(A — A) -1 , 
j = 1,2, are also integral operators of Hilbert-Schmidt class; in par- 
ticular, domi? D domA,, B is A,-compact [251 Ch. IV] and thus Aj- 
bounded with relative A,-bound [61 Lemma III. 2. 16]. As usual, I 
shall stand for the identity operator in L 2 (0, 1). 

Now the spectral problems (11. ip . (11. 2p and ( II. ip . ( II. 3p can be re- 
garded as the spectral problems for the quadratic operator pencils 
T 2 (p,r) and 7i(p, r) respectively, given by 

T j (p,r)(\):=\ 2 I-2\B-A j , j = 1,2, 

for A E C on the A-independent domain dom Tj(p,r) := domA,-. From 
the above-listed properties of the operators Aj and B we conclude that 
for every A E C, the operator Tj(p, r)(A) is well defined and closed on 
domTj(p, r). 

We now recall the following notions of the spectral theory of operator 
pencils, see [31]. 

4 



Definition 2.1. An operator pencil T is an operator- valued function 
on C. The spectrum cr(T) of an operator pencil T is the set of all A G C 
for which T(A) is not boundedly invertible, i.e. 

cr(T) = {A G C | G <r(T(A))}. 

A number A G C is called the eigenvalue of T if T(X)y = for some 
non-zero function y G domT, which is then the corresponding eigen- 
function. Finally, 

p(T):=C\a(T) 
is the resolvent set of an operator pencil T. 

It was shown in [38] that the spectra of the operator pencils Tj(p, r), 
j = 1,2, consist entirely of eigenvalues and that a(Tj(p, r)) are discrete 
subsets of C. There are examples of Tj(p,r) having non-real and/or 
non-simple eigenvalues; the latter means that the algebraic multiplicity 
of an eigenvalue can be greater than 1 [31]. In such a generic setting, 
the inverse spectral problem of reconstructing the potentials p and q 
in (II. ip is very complicated, and only partial existence results have been 
established so far. The approach we are going to use gives a complete 
solution in the special case where the spectra of Ti(p, r) and T 2 (p, r) 
are real and simple. A sufficient condition for this to hold is, e.g., that 
there is a /x* G R such that the operator Ti(p, r)(/i#) is negative [38] . 
Thus our standing assumption is that 

(A) the operator pencil Ti(p, r) is hyperbolic [3T], Ch. IV, §31], 
or, equivalently, there is a (i, G R such that the opera- 
tor Ti(p, r)(//#) is negative. 

Under assumption (A), the eigenvalues n G Z, of r) and A n , 
n G Z* := Z \ {0}, of T 2 (p, r) are all real, simple (see [38J ) and can be 
labelled in increasing order so that fi n and A n obey the asymptotics 

fl n = 7T (n - |) + p + fi n , \ n = 7in+p + A n , 

withpo : = Jo p( x ) dx and (/in), (A„) in £ 2 (see [37]) and almost interlace 
in the sense that < Afc < /i^+i when fcGZ* (see [13] and Section [3]). 
Thus the pair of spectra ((A n ), (/i n )) forms then an element of the 
set SD defined as follows. 

Definition 2.2. We denote by SD the family of all pairs (A, fi) of 
increasing sequences A := (A n ) ne ^* and /x := (/i n ) ng z of real numbers, 
which satisfy the following conditions: 

(i) asymptotics: there is an h G R such that 

(2.1) X n = 7m + h + X n , ^ n = 7r(n- |) +h + fi n , 

where (A n ) is a sequence in £ 2 (Z*) and (fl n ) is from £ 2 (Z); 
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(ii) almost interlacing: 
(2.2) /i k < \ k < /i fc+ i for every k eZ*. 

Take an arbitrary A G (/i ,yUi) and denote by A* the sequence A 
augmented with Ao- Then the almost interlacing condition (12. 2 p means 
that the sequences A* and li strictly interlace. 

Remark 2.3. Assume (A) and let the enumeration of the eigenvalues 
of Ti(p, r) and T 2 (p, r) agree with (12. 2ft ; then the number /z* in assump- 
tion (A) belongs to (/io,/Ui). Moreover, the operator 7\(p, r)(A) is neg- 
ative for every A G (^0,/Ui), so that (A) holds with every /i* G (/io,/ii). 
The same holds for T 2 (p, r)(/i); in particular, then T 2 (p, r)(A) is nega- 
tive for every A G (/i ,/ii). 

Denote by li and A the spectra of the operator pencils Ti(p, r) 
and T2(p,r) respectively. Then the inverse spectral problem of inter- 
est is to reconstruct the operator pencils Tj(p,r), j = 1,2, given the 
spectra li and A. Namely, we want to reconstruct the potential p and 
the (regularized) potential r (determining both the potential q = r' 
and the right-end boundary condition for the operator pencil 7i(p, r) 
in (11.31) ). Some properties of the sets li and A are given in |38j; our 
aim here is to give a complete characterization of these spectral data 
and to find an algorithm reconstructing the potentials p and r from the 
spectral data. 

It is easy to see that with p = 0, the spectral problems for the 
operator pencils Tj(0, r) become the usual spectral problems Ajy = \ 2 y 
for the Sturm-Liouville operators A,- with potential q in W r 2 " 1 (0, 1); 
see [IUjf¥T] . In this case parts of the sequences fx and A are redundant in 
the sense that /ii- n = — /i n and A_„ = — A n . In [TH] the authors proved 
that, just as in the regular case when q is integrable, the spectra (/^) nG pj 
of Ai and (A^) ne pj of A 2 uniquely determine the regularized potential r; 
see also [42J for an alternative treatment. These papers also suggest 
the algorithm of reconstructing the potential r from the spectra of A\ 
and A 2 . 

In the inverse spectral problem of interest we want to determine two 
real- valued potentials p and r of the operator pencils Tj(p, r), j = 1, 2. 
Since the information contained in the spectra of Tj(p,r) is twice as 
large as for the standard Sturm-Liouville operators one may hope that 
the inverse spectral problem of reconstructing p and r from the spectral 
data of Tj(p,r) is well posed. 

As was already mentioned, the spectra of the operator pencils Tj(p, r), 
j = 1,2, form an element of SD. Our main result states that, con- 
versely, every element of the set SD coincides with the spectral data 
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for some operator pencils Tj(p,r) under consideration and that these 
operator pencils are uniquely determined by their spectra. 

Theorem 2.4. Assume that a pair (A, n) of sequences of real numbers 
is an element of SD. Then there exist unique p,r G L2.r(0, 1) such 
that fi and A are the spectra of the pencils Ti(p, r) and T 2 (p, r) respec- 
tively. Moreover, the operator pencil Ti(p,r) satisfies assumption (A). 

The proof of this theorem is constructive and suggests the recon- 
struction algorithm determining the potentials p and r from the sets A 
and /x, see Section [6j 

As was already mentioned, our approach consists in reducing the 
spectral problems for operator pencils Tj(p, r), j = 1, 2, to the spectral 
problems for Dirac operators in L 2 (0, 1) X L 2 (0, 1) with appropriate 
boundary conditions, see Section [3J Using a suitable unitary gauge 
transformation we can reduce these Dirac operators to those in the 
"shifted" AKNS normal form. For AKNS Dirac operators, the direct 
and inverse spectral problems are well understood, see [3l[7|,|28]. Us- 
ing known methods we shall reconstruct the AKNS Dirac operators 
from the given spectra and then transform them to the Dirac opera- 
tors directly associated with the operator pencils of interest keeping 
the spectra unchanged. This will give the sought potentials p and r in 
an explicit form, see f)5.6p . 

3. Reduction to the Dirac system 

In this section we are going to show that under assumption (A) 
the spectral problems for the operator pencils Ti(p,r) and T 2 (p,r) 
can be reduced to the spectral problems for special Dirac operators. 
The reduction crucially relies on the fact that the equation y" = q*y 
with := q + 2p*p — pi and //* from (A) has a solution which is strictly 
positive on [0,1]. 

Observe first that this equation can be recast as £(y) = ply — 2p*py 
and treated as the first order linear system u[ = ru\ + u 2 , u' 2 = 
(2p*p — pi — r 2 )u\ — ru 2 , where U\ = y and u 2 = yW. Therefore 
for every complex a and b it possesses a unique solution satisfying the 
conditions y(l) = a and 2/^(1) = b. Let z denote the solution of the 
equation y" = q*y subject to the conditions 

2(1) = 1, Z [1] (1) = 0. 

We are going to show that z does not vanish on [0, 1]. 

Lemma 3.1. Under the standing assumption (A) the function z is 
strictly positive on [0,1]. 

7 



Proof. By assumption, the operator — Ti(p, r)(ju*) is uniformly posi- 
tive and thus such is the closure ti of its quadratic form. For y G 
dom Ti(p,r) we find by integration by parts that 

ti\y] = (-T 1 tp,r)( f ji,)y,y)= I \yf+ [ (2^p- ^)\y\ 2 -2 Re(ry', y). 

Jo Jo 

The quadratic form i [y] = \y'\ 2 considered on domT^, r) is clos- 
able and its closure to acts by the same formula on the domain 

dom t := {u G W 2 {0, 1) | w(0) = 0}. 

Since the quadratic form J^(2^p — n 2 )\y\ 2 — 2Re(ry\y) is relatively 
bounded with respect to the form io[y] with relative bound (cf. [15J ) , 
Theorem VI. 1.33 of [25] implies that the domain of ti coincides with 
dom t , i.e., 

domti = {ue W%(0, 1) | u(0) = 0}. 

Now assume that z has zeros on [0, 1] and denote the largest of them 
by Xq. Then the function 



111 x } 



0, X < Xq 

z(x), X > Xq 



belongs to the domain of ti and is not identically zero. It is straight- 
forward that ti[u] = 0, which contradicts positivity of — Ti(p, r) (//*). 
Therefore z does not vanish on [0, 1] yielding the assertion of the 
lemma. □ 

Setting v := z' / 'z, one easily verifies that g* = v' + v 2 ; hence g* is a 
Miura potential |24j and the differential expression —u" + q*u can be 
written in the factorized form, viz. 

(3.1) -„- + ,.„= -(^ + v){±-v)u. 

Since z satisfies (z' — rz)(l) = 0, we have that (v — r)(l) = 0. 

For A^/j t consider the functions u 2 '■— y and U\ := (y , —vy)/(X—fjL#). 
Then equation (ll.ip can be recast as the following first order system 
for u\ and w 2 , 

(3.2) u' 2 — vu 2 — (A — 

(3.3) -u[ - vui + 2pu 2 = (A + [i*)u 2 . 
Setting 



we see that the above system is the spectral problem £(P)u = Au for 
a Dirac differential expression £(P) acting in L 2 (0, 1) x L 2 (0, f ) via 

(3.5) £(P)u = J^ + Pu 

ax 

on suitable domain. 

Denote by @j(P), j = 1,2, the Dirac operators generated by £(P) 
on the domains 

dom^(P) := {u= {u 1 ,u 2 ) t G W^(0, 1) x W^O, 1) | u 2 (0) = ^(1) = 0}. 
It turns out that the spectra of ^j(P) and Tj(p, r) are closely related. 

Lemma 3.2. The spectra of the Dirac operators 2$j(P) and the oper- 
ator pencils Tj(p, r), j = 1,2, are related in the following way: 

(3.6) o-(® l (P))=o-(T 1 (p,r)), 

(3.7) a(^ 2 (P))=a(T 2 (p,r))U{^}. 

Moreover, y is an eigenfunction ofTj(p,r) corresponding to an eigen- 
value A if and only if u = (1/1,1*2)* with u\ := (y' — vy)/(\ — /i*) 
and U2 = y is an eigenfunction of @j (P) corresponding to A ^ //* . 

Proof. It is straightforward to see that A = /i* is an eigenvalue 
of the Dirac operator £^ 2 (P), the corresponding eigenfunction being 
u = ("Ui, , u 2 ) t with u\ = exp(— j v) and w 2 = 0. However, by Re- 
mark 12.31 under the standing assumption (A) the number A = /i* is 
not an eigenvalue of T 2 (p, r). The proof of (13.71) uses this fact and is 
analogous to the proof of the corresponding lemma in [T2|. The co- 
incidence of the spectra a(^i(P)) and a(Ti(p, r)) is proved along the 
same lines using the fact that the boundary condition = is 

equivalent to wi(l) = in view of the relation (v — r)(l) = 0. The last 
statement of the lemma is obtained by direct verification. □ 

4. Transformation operator 

Lemma 13.21 suggests that we can try to use the spectra [i and A of 
the operator pencils Ti(p,r) and T 2 (p,r) in order to find the related 
Dirac operators Stj(P), j = 1, 2. Having determined the potential P = 
(pij)fj =1 of £(P), we then identify the potentials p and r of the operator 

pencil Tj(p,r) as p := (p 22 +Pn)/2 and r := -p 12 - jl(p\ 2 -£>np 22 ). 

However, the classical inverse spectral theory reconstructs a Dirac 
operator with potential in the AKNS form or in other canonical form. 
Thus to find the Dirac operators directly associated with the opera- 
tor pencils Tj(p,r), i.e. with potentials of the form (13. 4p . we have to 
transform the obtained Dirac operators in canonical form keeping the 
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spectral data unchanged. This is done by means of the so-called trans- 
formation operators. In [14] such operators were constructed and some 
of their properties were discussed; see also [5j[28]. In this section we 
shall focus our attention on some further properties of the transforma- 
tion operators that are important for our reconstructing procedure. 

Assume that P and Q are 2x2 matrix- valued potentials in L 2 ((0, 1), A4 2 ) 
and set 

V := {(ui,^)* G ^(0,1) x 1^(0,1) | u 2 {0) = 0}. 

The transformation operator = 3£{P, Q) between the Dirac opera- 
tors £(P) and £(Q) on the set Do, i.e. the nontrivial operator satisfying 
the relation ^£{P)\i = £{Q)^u for all u e Vq, was constructed in [13] 
in the form 

(4.1) &u(x) = R(x)u(x) + / K(x,s)u(s)ds, 

Jo 

where R and K are 2x2 matrix-valued functions of one and two 
variables respectively. Under the normalization R(0) = I, the operator 
R is explicitly given j5] as 

(4 2) R(x) = e 8l[x) ( cos ° 2 ^ sin ^(a;) \ _ ^ «,,•,■)/-«,,.,■,./ 
^ ' ' ^ ' \ — sin6 l 2 (x) cos9 2 (x) 



with 



(4.3) 



9 1 (x) = ~J X tr[J(Q(s)-P(s))]d S , 

Ux) = \ f tr(Q(s) - P(s))ds. 
* Jo 



Existence of the transformation operator &(P, Q) is guaranteed by the 
following two theorems. 

Theorem 4.1 ([E]). Assume that P and Q are in L 2 ((0, 1), M 2 )- 
Then an operator 3£{P, Q) of the form ( 14. ip . with R obeying the con- 
dition R(0) = I and a summable kernel K, is a transformation operator 
for £(P) and £(Q) on the set T>q if and only if the matrix-valued func- 
tion R is given by (I4.2p - (l4.3p and the kernel K is a mild solution of 
the partial differential equation 

(4.4) Jd x K(x, y) + d y K{x, y)J = K(x, y)P(y) - Q(x)K(x, y) 

in the domain Q := {(x,y) |0<y<x<l} satisfying for < x < 1 
the boundary conditions 

(4.5) K(x, x)J — JK(x, x) = JR'(x) + Q(x)R(x) - R(x)P(x), 

(4.6) K 12 {x,0) = K 22 (x,0) = 0. 
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Theorem 4.2 ([2]). Assume that matrix-valued functions P andQ are 
in L 2 ((0, 1), A^)- Then the system f)4.4p -f l4TB|) has a unique solution in 
the sense of distributions; moreover, this solution belongs to L 2 (fl, M.2)- 

Throughout the rest of the paper, we shall assume that the matrix- 
valued potentials P and Q are Hermitian, i.e. that P*(x) = P(x) 
and Q*(x) = Q(x) a.e. on [0,1]. Then the corresponding Dirac op- 
erators £tj(P) and @j(Q), j = 1,2, are self-adjoint and have simple 
discrete spectra. We denote by A*(P) (resp. by A*(Q)) the spectrum 
of ^i(P) (resp. of 3!\{Q)) and by A(P) (resp. by A(Q)) the spectrum 
of ^(P) (resp. of (<?))• The sets i-i(P) and A(P) interlace and their 
elements can be labelled by n G Z so that /i n = tc (n — |) + | tr P + 

o(l) and A„ = irn + | J tr P + o(l) as |n| — > oo (see [28]); the spectra 
fj.(Q) and A(Q) have similar properties. 

Also, 3£ = <%~(P,Q) will stand for the transformation operator of 
the form (14. ip for the differential expressions i(P) and £(Q) on the 
domain V , with R given by (14. 2p and (14. 3p . We shall write = 
^ + J?T, where M\i(x) := R(x)u(x) is the operator of multiplication 
by R and 



is the corresponding integral operator. We observe that for Hermit- 
ian P and Q the functions %Q\ and 9 2 are real valued; in particular, the 
operator M is unitary. 

Remark 4.3. Let us note that if u G Vq, then the relation 



Theorem 4.4. Assume that the matrix potentials P and Q are Her- 
mitian. Then the spectra of the operators ^j(P) and S)j(Q), j = 1,2, 
coincide, i.e. l-i{P) = fi{Q) and A(P) = X(Q), if and only if the trans- 
formation operator S£(P, Q) for £(P) and £(Q) on the domain V only 
contains the unitary part M (i.e. = 0) and #2(1) = ?rn for some 
neZ. 

Proof. Sufficiency. By assumption, the transformation operator S£{P, Q) 
between £(P) and £{Q) only contains the multiplication operator 
Denote by Stj(P), j = 1, 2, the operators related to S'jiQ) by 




(£(P)-A)u = f 



holds if and only if for v := 3£vl and g := i?f f one gets 



(*(Q)-A)v = g. 



@j(P) = M^^j{Q)M, 3 = 1,2. 
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Then @j(P) are unitarily equivalent to @j(Q) and so a(@j(P)) = 
a(@j(Q)). By Remark H~3] @j(P) are Dirac operators acting via i^(P)u = 
£(P)u on the domains which consist of those u G L2((0,1),C 2 ) for 
which G dom^(Q). Recall that R(0) = I; also since # 2 (1) = im 
for some n G Z, we have that P(l) is a multiple of the identity 
matrix I. Therefore these u satisfy the same boundary conditions 
as ^u. Thus dom^(P) = dom^(Q) = dom^j(P). This means 
that %(P) = @j(P), and so fi(P) = n{Q) and A(P) = \{Q). 

Necessity. Firstly note that as the spectra of @j(P) and $>j(Q), 
j = 1,2, coincide, the asymptotics of the eigenvalues also coincide thus 
yielding 6 2 (1) = 7rn for some n G Z. To complete the proof it remains 
to show that Jff = 0. 

Recall that the norming constant corresponding to an eigenvalue A of 
the operator £^ 2 (P) is defined as ||u|| 2 , where u = (ui,U2) t is the eigen- 
function of @%(P) for A normalized by the initial conditions wi(0) = 1 
and «2(0) = [28J. It follows from [3] that the norming constants 
for S>2(P) are uniquely determined by the spectra A(P) and //(P) 
(although only AKNS potentials were treated in [3J, the arguments 
therein are valid for all self-adjoint Dirac operators). Thus from the 
statement of the theorem we conclude that the sets of norming con- 
stants for S^(P) and ^(Q) coincide. By Lemma 4.4. of [H] this gives 
that Jt = 0. □ 

Alternatively, to establish the necessity part of the above theorem, 
one can use arguments analogous to those in [5] , where a similar state- 
ment but under different assumptions on P and Q was proved. 

5. Reconstruction of the pencils 

In this section, we study the problem of reconstructing the poten- 
tials p and r of the corresponding operator pencils Tj(p,r), j = 1,2, 
from their spectra. Our aim is to prove Theorem 12 A\ i.e. that given an 
arbitrary element (A, n) of SD, there exist unique p and r such that \x 
is the spectrum of the operator pencil Ti(p, r) and A is that of T 2 (p, r). 

Fix therefore an arbitrary pair (A, /x) in SD. In what follows, A n , n G 
Z*, and fi n , n G Z, will stand respectively for elements of A and fi. Note 
that the enumeration of fi n and A n fixes the shift h in their asymp- 
totics (12.11) . We set /1* to be the middle point of the interval (/iQ,/ii), 
i.e. //* := (/io + £*i)/2, and augment the sequence (X n )nez* with an 
element Ao = //*; denote obtained sequence by A*. 

Recall now some facts from the inverse spectral theory for Dirac 
operators, which we shall use in our procedure, see [3JC3I28]. Consider 
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the set Qq of 2 x 2 matrix-valued functions of the AKNS normal form, 
namely 



It is known that the operators @i(Qo) and ^(Qo) with potential Q 
from Q are self-adjoint and have simple discrete spectra. Their eigen- 
values can be enumerated in increasing order as X n (Qo) and fi n (Qo), 
n G Z, respectively so that they satisfy the interlacing condition 

Hn{Qo) < K(Qo) < AWl(Qo) 
and obey the asymptotics 



It is also known (see j3]) that for two sequences (/x n ) and (A n ) having 
the above properties there exists a unique potential Q from Q such 
that the spectrum of @i(Qo) is (/x n ) and that of & 2 {Qo) is (A n ). 

The above sequences /x and A* differ from the spectra for Dirac 
operators with potentials of normal AKNS form only by the shift h in 
their asymptotics, see (12.11) . For h G R, we denote by 



the set of h-shifted AKNS matrix potentials; then the following result 
holds true. 

Proposition 5.1. For an arbitrary pair (A, /x) from SD take /x* : = 
(/io+/ii)/2 and denote by A* the augmentation of X with \ = /x*. Then 
there exists a unique potential Q G Qh such that /x is the spectrum of 
the operator Q>\{Q) and A* is that of @ 2 (Q). 

For an arbitrary P G L 2 ((Q, 1), M.2), we introduce the set 



Iso(P) := {P G L 2 ((0,1),M 2 ) | H(P) = fi(P), A(P) = X(P)} 



of isospectral potentials and for /x G R denote by P M the set of all 
potentials of the form (13.41) . i.e. 

'Pfj, ■= {P = (Pij)lj=l I Pij € L 2 ,r(0, 1), Pli = /X, Pl2 = P2l}- 

Now we are going to prove that for every /x G R and Q £ Qh there ex- 
ists a unique potential P from lso(Q) fTP M . We start with the following 



(5.1) 




An(Qo) =7™ + An(Qo), 
^n(Qo) = 7T (n - |) + /X„(Q ) 

with £ 2 (Z)-sequences (\ n (Qo)) and (/x n (Qo))- 



:= {Qo + /*/ I Q e So} 



lemma. 
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Lemma 5.2. Suppose that Q G Qh, i-e. 

q = ( Qi + h q 2 

Then for every fixed /i 6 R there exists a unique potential P G such 
that, with 1Z being the operator of multiplication by the matrix-valued 
function R of (14. 2 j) - (14.31) . one /ias 

(5.2) ^£(P) = £(Q)^. 

Proof. First we note that the matrix R given by (14.21) commutes with J, 
whence 

r~ 1 (j-^- + q)r = j-^- + r~ 1 jr' + r~ 1 qr. 

\ ax / ax 
Thus the relation (15. 2p requires that 
(5.3) 

P = R^JR' + R~ l QR 

- (h - fl'M 4- 0' 7- 4- ( 9i cos 20 2 - g 2 sin 20 2 <?i sin 20 2 + q 2 cos 20 2 \ 
1 2 > + 1 + ^ gi sin 20 2 + 92 cos 20a ~9i cos20 2 + q 2 sin20 2 / 

Now we fix an arbitrary // G R and observe that the potential P of (15.31) 
belongs to if and only if 9\{x) is identically zero and the following 
equality holds: 

(5.4) -#2 + 9i cos 20 2 - q 2 sin 20 2 + h = //. 

Clearly, there exists a unique solution of the above equation satisfying 
the initial condition 

(5.5) 2 (O) = 0. 

This solution and 6\ = verify the relation (14.31) with the given Q and 
with P of (15.31) . The potential P is explicitly given by (15. 3p with Q\ = 
and 2 as above and by construction it belongs to P M and satisfies (15. 2p . 
The proof is complete. □ 

The next result states existence and uniqueness of the potential P 
from Vfj, belonging to Iso(<5)- 

Theorem 5.3. Let ^, h G R and assume that Q G Qh is such that 
X — ii is an eigenvalue of the Dirac operator @ 2 (Q). Then there exists 
a unique P G P M belonging to Iso(Q), i.e. Iso(Q) H P M = {P}. 

Proof. Observe first that if P G Iso(Q), then by Theorem 14.41 the trans- 
formation operator between £(P) and £(Q) on T>q is just the operator & 
of multiplication by the matrix- valued function R given by (I4.2l) - (l4.3p . 
If, moreover, P should belong to P M , then by Lemma l5~2l P should be 
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v (l) = i2(l)uo(l) 



given by (15.31) with &x = and 82 solving fl5.4p -f l5"3j) . It remains to 
prove that the potential given by fl5.3p - fl5.5p indeed belongs to Iso(Q). 
By Theorem I4.4[ to this end it suffices to show that #2(1) — for 
some fiGZ. 

Consider the operator %{P) ■= M- l 9 2 {Q)M. Clearly, 9 2 {P) is 
a Dirac operator acting via ^(-P)u = £(P)u on the domain con- 
sisting of those u for which TZu belongs to dom^ 2 (Q)- By con- 
struction, a[S>i{Py) = a(^2(Q)) and thus A = ji is in the spectrum 
of 3>2(P)- Denote by u = (141,1X2)* the eigenfunction of i^ 2 (P) corre- 
sponding to A = fi. Then v = (fi,^)* := ^u is an eigenfunction 
of & 2 {Q) corresponding to A, so that t>2(0) = t>2(l) = 0. The equal- 
ity £(P)u = fiu yields the relation u' 2 — vu 2 = 0. Since R(0) = I, we 
have 1*2(0) = ^2(0) = 0, which together with the above relation gives 
that U2 is identically zero. Therefore 

ui(l)cos0 2 (l) 
- Ul (l)8in0 2 (l) 

But i> 2 (l) = and «i(l) 7^ giving that 02(1) = ^n, n E Z. □ 

With these preliminaries at hand we can prove the main results of 
the paper. 

Proof of Theorem \2.4\ Existence. Given a pair (A, fi) from SD, we 
set /i* := (fiQ + fii)/2 and augment A with A := //* to obtain A*. 
Using Proposition 15. 11 we find the potential Q in Q h such that fx is the 
spectrum of the operator @\{Q) and A* is that of ^(Q)- Theorem 15.31 
gives that for this Q we can find P G belonging to Iso(Q), i.e. such 
that li and A* are respectively the spectra of @i(P) and @ 2 (P)- In 
view of Lemma [3721 this means that ll and A are the spectra of Ti(p, r) 
and T 2 (p, r) with 

fx a\ P22+P11 f \ 2 \ 

(5-6) p = , r := -p 12 - J (p 12 - p n p 22 ). 

The almost interlacing of A and fi gives that T\(p, r) is hyperbolic 
(see |38j), i.e. it satisfies assumption (A). 
This completes the proof of existence. 

Uniqueness. Suppose there are two pairs of potentials p, r and p, f 
such that a(Ti(p, r)) = <r(Ti(p, f)) =: li and cr(T 2 (p, r)) = a{T2{p 1 r)) =: 
A. Then (see [38]) for every p e (// , /^i) the operators Ti(p, r)(p) 
and Ti(p, r)(p) are negative, i.e. the assumption (A) holds for operator 
pencils 7\(p, r) and Tx(p, f). 

Set p* := (p + A*i)/2; then the pencils Tj(p,r) and Tj(p,r), j = 
1,2, lead to two Dirac operators with potentials P and P in T 7 ^ as 
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explained in Section [3] and to two Dirac operators with potentials Q 
and Q in Qh, with the same h as explained at the beginning of this 
section. By construction, fi(Q) = fJ>(Q) = A* and \(Q) = \(Q) = 
A U {/i*}, whence Q = Q by Proposition 15.11 By Lemma [3.21 we have 
that p(P) = fj,(P) and A(P) = A(P) and thus P and P belong to 
the same isospectral set Iso(Q). In view of Theorem 15.31 this yields 
that P = P. Therefore, p = p and f = f, which completes the proof. 

□ 

6. Reconstructing algorithm 

To summarize let us formulate the algorithm of reconstruction of the 
operator pencils Ti(p, r) and T 2 (p,r). 

Suppose we have a pair of sequences (A, n) from SD. Then the 
reconstructing procedure consists of the following steps: 

(1) Take /x* := (/ii + /i2)/2 and augment the given sequence A with 
an element Ao = /U*; 

(2) Use the obtained spectral data to construct the Dirac operators 
V X {Q) and V 2 (Q) with the potential Q of the /i-shifted AKNS 
form; 

(3) Find the corresponding potential P G fl Iso(<5) using (15.31) 
and (J53D; 

(4) Compute the potentials p and r using formulae (15.61) . 

We conclude this paper with some remarks. Firstly, note that the de- 
scribed approach can be used to reconstruct energy-dependent Sturm- 
Liouville equations subject to other boundary conditions or under dif- 
ferent smoothness assumptions on the potentials, e.g. in the case when p 
and r are from L s (0, 1) with s > 1 or belong to W^|(0, 1) with s > 0, 

cf. [tang. 

Secondly, the considered method is also applicable to the inverse 
problems of reconstructing the quadratic pencil from different sets 
of spectral data, for example, the spectrum and the norming con- 
stants [Ej or Hochstadt-Lieberman mixed data (cf. [12]). Finally, one 
can establish an analogue of the Hochstadt-type perturbation formu- 
las when finitely many eigenvalues in one spectrum have been changed 

(cf- HU). 
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